Next-generation sequencing identified that RET variation associates with lymph node metastasis and the immune microenvironment in thyroid papillary carcinoma.

BACKGROUND: To date, although most thyroid carcinoma (THCA) achieves an excellent prognosis, some patients experience a rapid progression episode, even with differentiated THCA. Nodal metastasis is an unfavorable predictor. Exploring the underlying mechanism may bring a deep insight into THCA. METHODS: A total of 108 THCA from Chinese patients with next-generation sequencing (NGS) were recruited. It was used to explore the gene alteration spectrum of THCA and identify gene alterations related to nodal metastasis in papillary thyroid carcinoma (PTC). The Cancer Genome Atlas THCA cohort was further studied to elucidate the relationship between specific gene alterations and tumor microenvironment. A pathway enrichment analysis was used to explore the underlying mechanism. RESULTS: Gene alteration was frequent in THCA. BRAF, RET, POLE, ATM, and BRCA1 were the five most common altered genes. RET variation was positively related to nodal metastasis in PTC. RET variation is associated with immune cell infiltration levels, including CD8 naïve, CD4 T and CD8 T cells, etc. Moreover, Step 3 and Step 4 of the cancer immunity cycle (CIC) were activated, whereas Step 6 was suppressed in PTC with RET variation. A pathway enrichment analysis showed that RET variation was associated with several immune-related pathways. CONCLUSION: RET variation is positively related to nodal metastasis in Chinese PTC, and anti-tumor immune response may play a role in nodal metastasis triggered by RET variation.

VDR is a potential prognostic biomarker and positively correlated with immune infiltration: a comprehensive pan-cancer analysis with experimental verification.

The vitamin D receptor (VDR) is a transcription factor that mediates a variety of biological functions of 1,25-dihydroxyvitamin D3. Although there is growing evidence of cytological and animal studies supporting the suppressive role of VDR in cancers, the conclusion is still controversial in human cancers and no systematic pan-cancer analysis of VDR is available. We explored the relationships between VDR expression and prognosis, immune infiltration, tumor microenvironment, or gene set enrichment analysis (GSEA) in 33 types of human cancers based on multiple public databases and R software. Meanwhile, the expression and role of VDR were experimentally validated in papillary thyroid cancer (PTC). VDR expression decreased in 8 types and increased in 12 types of cancer compared with normal tissues. Increased expression of VDR was associated with either good or poor prognosis in 13 cancer types. VDR expression was positively correlated with the infiltration of cancer-associated fibroblasts, macrophages, or neutrophils in 20, 12, and 10 cancer types respectively and this correlation was experimentally validated in PTC. Increased VDR expression was associated with increased percentage of stromal or immune components in tumor microenvironment (TME) in 24 cancer types. VDR positively and negatively correlated genes were enriched in immune cell function and energy metabolism pathways, respectively, in the top 9 highly lethal tumors. Additionally, VDR expression was increased in PTC and inhibited cell proliferation and migration. In conclusion, VDR is a potential prognostic biomarker and positively correlated with immune infiltration as well as stromal or immune components in TME in multiple human cancers.

A distinct tumor microenvironment makes anaplastic thyroid cancer more lethal but immunotherapy sensitive than papillary thyroid cancer.

Both anaplastic thyroid cancer (ATC) and papillary thyroid cancer (PTC) originate from thyroid follicular epithelial cells, but ATC has a significantly worse prognosis and shows resistance to conventional therapies. However, clinical trials found that immunotherapy works better in ATC than late-stage PTC. Here, we used single-cell RNA sequencing (scRNA-Seq) to generate a single-cell atlas of thyroid cancer. Differences in ATC and PTC tumor microenvironment components (including malignant cells, stromal cells, and immune cells) leading to the polarized prognoses were identified. Intriguingly, we found that CXCL13+ T lymphocytes were enriched in ATC samples and might promote the development of early tertiary lymphoid structure (TLS). Last, murine experiments and scRNA-Seq analysis of a treated patient's tumor demonstrated that famitinib plus anti-PD-1 antibody could advance TLS in thyroid cancer. We displayed the cellular landscape of ATC and PTC, finding that CXCL13+ T cells and early TLS might make ATC more sensitive to immunotherapy.

Single-cell RNA-seq reveals MIF-(CD74 + CXCR4) dependent inhibition of macrophages in metastatic papillary thyroid carcinoma.

BACKGROUND: The mechanism promoting papillary thyroid carcinoma (PTC) metastasis remains unclear. We aimed to investigate the potential metastatic mechanisms at a single-cell resolution. METHODS: We performed single-cell RNA-seq (scRNA-seq) profiling of thyroid tumour (TT), adjacent normal thyroid (NT) and lymph node metastasized tumour (LN) from a young female with PTC. Validation of our results was conducted in 31 tumours with metastasis and 30 without metastasis. RESULTS: ScRNA-seq analysis generated data on 38,215 genes and 0.14 billion transcripts from 28,839 cells, classified into 18 clusters, each annotated to represent 10 cell types. PTC cells were found to originate from epithelial cells. Epithelial cells and macrophages emerged as the strongest signal emitters and receivers, respectively. After reclustering epithelial cells and macrophages, our analysis, incorporating gene set variation analysis (GSVA), SCENIC analysis, and pseudotime trajectory analysis, indicated that subcluster 0 of epithelial cells (EP_0) showed a more malignant phenotype, and subclusters 3 and 4 of macrophages (M_3 and M_4) demonstrated heightened activity. Further analysis suggested that EP_0 may suppress the activity of M_3 and M_4 via MIF - (CD74 + CXCR4) in the MIF pathway. After analysing the expression of the 4 genes in the MIF pathway in both the TCGA cohort and our cohort (n = 61), CD74 was identified as significantly overexpressed in PTC tumours particularly those with lymph node metastasis. CONCLUSION: Our study revealed that PTC may facilitate lymph node metastasis by inhibiting macrophages via MIF signalling. It is suggested that malignant PTC cells may suppress the immune activity of macrophages by consistently releasing signals to them via MIF-(CD74 + CXCR4).

BRAFV600E and BRAF-WT Specific Antitumor Immunity in Papillary Thyroid Cancer.

One feature of papillary thyroid cancer (PTC) is the frequently present somatic BRAFV600E mutation. PTCs are also characterized by a lymphocytic infiltration, which may correlate with an improved clinical outcome. The objective of the study was the characterization of BRAFV600E specific anti-immunity in PTC patients and correlation analyses with the clinical outcome. Fourteen HLA A2 positive PTC patients were included into the study of whom tumor tissue samples were also available. Of those, 8 PTC patients revealed a somatic BRAFV600E mutation. All PTC patients were also MHC class II typed. Tetramer analyses for detection of MHC class I and MHC class II-restricted, BRAFV600E epitope-specific T cells using unstimulated and peptide-stimulated T cells were performed; correlation analyses between MHC phenotypes, T cell immunity, and the clinical course were performed. In regard to unstimulated T cells, a significantly higher amount of BRAFV600E epitope specific T cells was detected compared to a control tetramer. Importantly, after overnight peptide stimulation a significantly higher number of BRAFV600E positive and BRAF WT epitope-specific T cells could be seen. In regard to the clinical course, however, no significant differences were seen, neither in the context of the initial tumor size, nor in the context of lymph node metastases or peripheral metastastic spread. In conclusion, we clearly demonstrated a BRAF-specific tumor immunity in PTC-patients which is, however, independent of a BRAFV600E status of the PTC patients.

Low expression of TFF3 in papillary thyroid carcinoma may correlate with poor prognosis but high immune cell infiltration.

Background: Papillary thyroid carcinoma (PTC) is one of the most common endocrine malignancies and has a favorable prognosis. However, optimal treatments and prognostic markers have not been clearly identified. Methods: Gene expression data from primary PTC were downloaded from the Gene Expression Omnibus database and subjected to two analyses of differentially expressed genes (DEGs), followed by intersecting individual and integrated DEGs analyses as well as gene set enrichment analysis. Analysis of data from Sequence Read Archive and The Cancer Genome Atlas, immunohistochemistry and qRT-PCR of TFF3 were performed to validate the results. Finally, the relationship between gene expression and disease-free survival as well as immune cell infiltration were investigated. Results: Six critical DEGs and several tumor-enriched signaling pathways were identified. Immunohistochemistry and qRT-PCR validated the low expression of TFF3 in PTC. TFF3 and FCGBP are coexpressed in PTC, and patients with lower gene expression had worse disease-free survival but higher immune cell infiltration. Conclusion: TFF3 was significantly underexpressed and may function with FCGBP synergistically in PTC.

Lay abstract Thyroid cancers are some of the most common endocrine malignancies. However, the optimal treatments and prognostic markers have not been clearly identified. We identified six critical differentially expressed genes and several tumor-enriched signaling pathways in papillary thyroid carcinoma, and found that TFF3 was the most underexpressed gene, as validated by experiment. In addition, TFF3 and FCGBP worked synergistically and may mark prognosis and tumor immune cell infiltration, which may benefit patients with papillary thyroid carcinoma by providing early indication and prompting further basic investigation.

Proteomics and Organoid Culture Reveal the Underlying Pathogenesis of Hashimoto's Thyroiditis.

Hashimoto's thyroiditis (HT) is an autoimmune disease, and its incidence continues to rise. Although scientists have studied this disease for many years and discovered the potential effects of various proteins in it, the specific pathogenesis is still not fully comprehended. To understand HT and translate this knowledge to clinical applications, we took the mass spectrometric analysis on thyroid tissue fine-needle puncture from HT patients and healthy people in an attempt to make a further understanding of the pathogenesis of HT. A total of 44 proteins with differential expression were identified in HT patients, and these proteins play vital roles in cell adhesion, cell metabolism, and thyroxine synthesis. Combining patient clinical trial sample information, we further compared the transient changes of gene expression regulation in HT and papillary thyroid carcinoma (PTC) samples. More importantly, we developed patient-derived HT and PTC organoids as a promising new preclinical model to verify these potential markers. Our data revealed a marked characteristic of HT organoid in upregulating chemokines that include C-C motif chemokine ligand (CCL) 2 and CCL3, which play a key role in the pathogenesis of HT. Overall, our research has enriched everyone's understanding of the pathogenesis of HT and provides a certain reference for the treatment of the disease.

ANXA1 (Annexin A1) regulated by MYC (MYC proto-oncogene) promotes the growth of papillary thyroid carcinoma.

Thyroid carcinoma is one of the most common endocrine malignancies, in which papillary thyroid carcinoma (PTC) is the main pathotype. ANXA1 plays a significant role in many cancer types, but how it works in PTC has not been identified. MYC is a common transcript factor involved in tumorigenesis, development, invasion, and metastasis. The relation between ANXA1 and MYC has not been proved in PTC. In this study, firstly, we analyzed the expression and prognostic value of ANXA1 in pan-cancer using the data from the UCSC database. Then we explore the role of ANXA1 in PTC, including expression, prognostic value, and immune infiltration. In addition, we evaluated the relation between ANXA1 and the transcription factor MYC. Finally, we identified the expression of ANXA1 and MYC and then evaluated their function associated with proliferation and apoptosis in PTC cell lines by CCK8 proliferation and flow cytometry apoptosis experiment. We found that ANXA1 is up-regulated in PTC comparing with normal patients. High expression of ANXA1 was associated with adverse overall survival of PTC. ANXA1 may be regulated by MYC to promote the proliferation of PTC. MYC may regulate the expression of ANXA and thus affect the proliferation of PTC.

Development and validation of a ferroptosis-related prognostic model for the prediction of progression-free survival and immune microenvironment in patients with papillary thyroid carcinoma.

BACKGROUND: Ferroptosis is an iron-dependent and regulated cell death that has been widely reported in a variety of malignancies. The overall survival of papillary thyroid cancer (PTC) is excellent, but the identification of patients with poor prognosis still faces challenges. Nevertheless, whether ferroptosis-related genes (FRGs) can be used to screen high-risk patients is not clear. METHODS: We obtained the clinical data of patients with PTC and FRGs from the UCSC Xena platform and the FerrDb respectively. Differentially expressed genes (DEGs) of FRGs were obtained from the entire The Cancer Genome Atlas (TCGA). Subsequently, the entire TCGA dataset was randomly split into two subsets: training and test datasets. Based on DEGs, we constructed a predictive model which was tested in the test dataset and the entire TCGA dataset to predict progression-free survival (PFS). Patients were categorized into high- or low-risk groups based on their median risk score. We analyzed differences in some aspects, including pathway enrichment analysis, single-sample Gene Set Enrichment Analysis (ssGSEA), tumor microenvironment (TME), human leukocyte antigen (HLA) genes, and tumor mutation burden (TMB) analyses, between high-risk and low-risk groups. RESULTS: A predictive model with three FRGs (HSPA5, AURKA, and TSC22D3) was constructed. Patients in the high-risk group had worse PFS compared with patients in the low-risk group. Functional analysis results revealed that ssGSEA, immune cell infiltration, TME, HLA, and TMB were closely associated with ferroptosis. CONCLUSION: The prognostic model constructed in this study can effectively predict PFS for patients with PTC.

T cell exhaustion is associated with the risk of papillary thyroid carcinoma and can be a predictive and sensitive biomarker for diagnosis.

BACKGROUND: The incidence of papillary thyroid carcinoma (PTC) has been steadily increasing over the past decades. Hashimoto's thyroiditis (HT) is the most common autoimmune disease, and is related to the pathogenesis of PTC. Programmed death-1 (PD-1) is currently used for the treatment of PTC, but there are very few studies on the clinical value of PD-1 in the diagnosis and targeted therapy of PTC. METHODS: The expression of T, B, NK cells and PD-1 in the peripheral blood of 132 patients with PTC (PTC group), 48 patients with nodular goiter (NG group) and 63 healthy subjects (HP group) were detected by flow cytometry. The expression of plasma T3, T4, FT3, FT4, TSH, TGAb and TPO was detected by chemiluminescence immunoassay. Among 132 PTC, 49 PTC&HT and 83 PTC&noHT were included. Among 48 NG, 10 NG&HT and 38 NG&noHT were included. The expressions of programmed death- ligand1(PD-L1) in tumor tissues of PTC group and thyroid tissues of NG group, PD-1 and CD3 in tumor infiltration lymphocyte (TIL) were detected by immunohistochemistry. RESULTS: The expression of FT3, TGAb, CD3+PD-1+, CD3+CD4+PD-1+ and CD3+CD8+PD-1+ in PTC and NG was significantly higher than that in the HP group. Moreover, CD3+PD-1+, CD3+CD4+PD-1+ and CD3+CD8+PD-1+ expression had significant differences between the PTC group and the NG group. In addition, the expression of TGAb, TPO, CD3+PD-1+, CD3+CD4+PD-1+ and CD3+CD8+PD-1+ in PTC&HT group was significantly higher than that in the PTC&noHT group. While, the expression of B cells, CD3+PD-1+, CD3+CD4+PD-1+ and CD3+CD8+PD-1+ in PTC&HT group was higher than that in NG&HT group. PD-1 showed a significant correlation with PTC lymph node metastasis. CD3+PD-1+ and CD3+CD4+PD-1+ was higher in N1 stage than in N0 stage. Immunohistochemical results showed that the expression of PD-1, CD3 and PD-L1 in PTC was significantly higher than that in NG. CONCLUSIONS: T cell exhaustion might act as a biomarker for the differential diagnosis of PTC and NG. Patients with PTC&HT have obvious T cell exhaustion and increased expression of PD-1, PD-L1.Targeting the PD-1/PD-L1 pathway could be a new approach to prevent malignant transformation from HT to PTC&HT in the future.

Prevalence of Hashimoto Thyroiditis in Adults With Papillary Thyroid Cancer and Its Association With Cancer Recurrence and Outcomes.

Importance: Hashimoto thyroiditis (HT) has been suggested to be associated with papillary thyroid cancer (PTC) development. However, its association with PTC progression remains unclear. Objective: To examine the association between HT and PTC presentation and outcomes. Design, Setting, and Participants: This retrospective cohort study included a review of patients aged 18 to 75 years who had pathologically confirmed PTC treated at a single center in China from January 1, 2001, to December 31, 2014. Data analysis was performed from November 1 to December 31, 2020. Exposures: Coexistent HT was defined according to evaluation of postoperative paraffin sections. Main Outcomes and Measures: The primary outcome was the association of HT with PTC-related mortality, assessed using Cox proportional hazards regression models. The secondary outcome was the association of HT with aggressive characteristics and structural recurrence of PTC, assessed using logistic regression and Cox proportional hazards regression with and without adjustment for related factors. Results: Of 9210 patients with PTC (mean [SD] age, 43.6 [12.0] years; 6872 [75%] women) included in the analysis, 1751 (19%) had HT. In the logistic regression model, HT was negatively associated with frequencies of primary tumor size of 4 cm or greater (adjusted odds ratio [aOR], 0.20; 95% CI, 0.12-0.33; P < .001), gross extrathyroidal extension (aOR, 0.44; 95% CI, 0.36-0.54; P < .001), extranodal extension (aOR, 0.66; 95% CI, 0.55-0.80; P < .001), and distant metastasis (aOR, 0.17; 95% CI, 0.04-0.71; P = .02). After a median follow-up of 85 months (range, 12-144 months), 131 PTC-related deaths were identified in the cohort; 2 patients who died had HT. Patients with HT had significantly superior outcomes compared with patients without HT in terms of unadjusted 10-year disease-specific survival (99.9% vs 96.6%; log-rank P < .001) and recurrence-free survival (92.0% vs 87.6%; log-rank P = .001). After adjusting for sex, age, primary tumor size, extrathyroidal extension, lymph node metastasis, distant metastasis, extent of surgery, and radioactive iodine ablation, HT was associated with decreased PTC-related mortality (hazard ratio [HR], 0.19; 95% CI, 0.05-0.76; P = .02). Stratified analysis showed that HT was associated with less frequent structural recurrence in patients with extrathyroidal extension (HR, 0.52; 95% CI, 0.38-0.71; P < .001; P = .002 for interaction) or after total thyroidectomy (HR, 0.50; 95% CI, 0.35-0.69; P < .001; P = .009 for interaction). Conclusions and Relevance: In this cohort study, patients with coexistent HT had less aggressive characteristics at presentation and better outcomes of PTC than did patients without HT. The findings suggest that autoimmune thyroiditis has a protective role in association with thyroid cancer.

Tumor Mutation Burden Predicts Relapse in Papillary Thyroid Carcinoma With Changes in Genes and Immune Microenvironment.

Background: The risk factors of papillary thyroid carcinoma (PTC) recurrence are meaningful for patients and clinicians. Tumor mutation burden (TMB) has been a biomarker for the effectiveness of immune checkpoint inhibitor (ICI) and prognosis in cancer. However, the role of TMB and its latent significance with immune cell infiltration in PTC are still unclear. Herein, we aimed to explore the effect of TMB on PTC prognosis. Material and Methods: RNA-seq and DNA-seq datasets of PTC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Gene Ontology (GO) and gene set enrichment analysis (GSEA 4.0.1) were applied further to explore potential differences in PTC patients' biological functions. The differentially expressed genes (DEGs) and immune microenvironment between the high and low TMB groups were determined. Results: TMB had the highest AUC score than other clinical indicators in ROC analysis on recurrence-free survival, and a higher TMB score was related to a worse prognosis. Further, GSEA showed a higher level of oxidative phosphorylation (OXPHOS) in the high TMB group, and four genes correlated with recurrence-free survival rate were identified. The abundance of CD8+ T cells and M1 macrophages in the high TMB group was significantly lower than that in the low TMB group. Conclusions: Our study found that TMB was a better predictor variable at evaluating the risk of PTC recurrence. Moreover, TMB-related genes conferred dramatically correlated prognosis, which was worth exploring in guiding postoperative follow-up and predicting recurrence for PTC patients.

Diagnostic and prognostic value of tumor-infiltrating B cells in lymph node metastases of papillary thyroid carcinoma.

Lymph node metastases are strongly associated with unfavorable prognosis in papillary thyroid carcinoma (PTC) patients. However, there are few sensitive or specific indicators that can diagnose or predict lymph node metastases in PTC. The objective of our study was to identify reliable indicators for the diagnosis and prediction of lymph node metastases of PTC. The PTC data set was obtained from The Cancer Genome Atlas (TCGA) cohort. Information on tumor-infiltrating immune cells in PTC was acquired using single-sample gene set enrichment analysis (ssGSEA). Then, the progression-free survival (PFS) rates of PTC patients were evaluated by Kaplan-Meier curves. A tissue microarray including 58 normal thyroid tissues and 57 PTC tissues was processed for CD19 immunohistochemistry staining. Finally, evaluation of phenotype permutations was performed using gene set enrichment analysis (GSEA). There was an appreciable association between immune infiltration and lymph node metastases in PTC. Among those immune cells, B cells and cytotoxic cells showed significant predictive accuracy for lymph node metastases in PTC. Tumor-infiltrating B cells and NK cells were associated with favorable prognosis, while tumor-associated NK CD56bright cells were correlated with poor prognosis in PTC patients. IHC analyses of PTC further confirmed a notably negative correlation between B cell infiltration and lymph node metastases in PTC. Additionally, mutations in BRAF, a dominant cause of tumor mutation burden (TMB), were positively correlated with reduced B cell infiltration and lymph node metastases in PTC. GSEA revealed that epithelial-mesenchymal transition, IL-6/JAK/STAT3 signaling, the inflammatory response, and TNF-α signaling via the NFκB pathway were remarkably suppressed pathways in patients with BRAF mutations. Tumor-associated lymphocytic infiltration, especially B cell infiltration, provides diagnostic and prognostic value for lymph node metastases in PTC.

T Cell Receptor Beta-Chain Profiling of Tumor Tissue, Peripheral Blood and Regional Lymph Nodes From Patients With Papillary Thyroid Carcinoma.

Objective: To study the characteristics of the T cell receptor (TCR) repertoire in cancer tissue, peripheral blood and regional lymph nodes (LNs) from patients with papillary thyroid carcinoma (PTC). Methods: PTC tissue, peripheral blood mononuclear cells (PBMCs) and regional LNs of six patients with papillary thyroid carcinoma were harvested. T cell receptor beta-chain (TCRß) profiling was performed though high-throughput sequencing (HTS), and IMonitor, MiXCR and VDJtools were used to analyze the characteristics of the TCR repertoire. Results: The results of IMonitor and those of MiXCR and VDJtools were very similar. The unique CDR3 of TCRß from LNs was higher than that of PBMCs, and the CDR3 of TCRß from LNs was higher than that of PTC tissue. Shannon's diversity index, D50, inverse Simpson index_mean and normalized Shannon's diversity index_mean of CDR3 from LNs were higher than those of PTCs and PBMCs. The HEC (high expansion clones) rate of CDR3 sequences at the amino acid level in PTC tissue was higher than that of PBMCs, which was higher than that of LNs. The V-J HEC rate of CDR3 was highest in PTC tissue, followed by PBMCs and LNs. Conclusion: TCR CDR3 profiling showed differences among and within the PBMCs, PTC tissues and regional LNs of PTC, including unique CDR3, CDR3 HEC at the amino acid level, CDR3 V-J HEC at the amino acid level, Shannon's diversity index and D50. The TCRß repertoire of PTC tissue, peripheral blood and regional LNs of PTC provide a reference for further study of immunity mechanisms against PTC.

Development and validation of a three-immune-related gene signature prognostic risk model in papillary thyroid carcinoma.

PURPOSE: Increasing evidence indicates that there is a correlation between papillary thyroid carcinoma (PTC) prognosis and the immune signature. Our goal was to construct a new prognostic tool based on immune genes to achieve more accurate prognosis predictions and earlier diagnoses of PTC. METHODS: The 493 PTCs samples and 58 tumor-adjacent normal tissues were obtained from The Cancer Genome Atlas database (TCGA). Immune genes were obtained from the ImmPort database. First, this cohort was randomly divided into training cohort and testing cohort. Second, the differentially expressed (DE) immune genes from the training set were used to construct the prognostic model. Then, the testing and entire data cohorts were used to validate the model, and the data were analyzed to determine the correlation of the clinical prognostic model with immune cell infiltration and expression profiles of human leukocyte antigen (HLA) genes. Finally, an analysis of the gene ontology (GO) annotation was performed. RESULTS: A total of 189 upregulated and 128 downregulated DE immune genes were identified. We developed and validated a three-immune gene model for PTC that includes Hsp70, NOX5, and FGF23. This model was demonstrated to be an independent prognostic variable. In addition, the overall immune activity of the high-risk group was higher than that of the low-risk group. CONCLUSIONS: We developed and validated a three-immune gene model for PTC that includes HSPA1A, NOX5, and FGF23. This model can be used as a validated tool to predict outcomes in PTC.

Potential role of M2 TAMs around lymphatic vessels during lymphatic invasion in papillary thyroid carcinoma.

The aim of this study was to examine whether lymphatic invasion in papillary thyroid carcinoma (PTC) occurs when tumour-associated macrophages (TAMs) injure lymphatic vessels together with cancer cells. While there was no difference in the lymphatic vessel density in PTC and follicular thyroid carcinoma (FTC), the number of TAMs around the lymphatic vessels was increased in PTC compared to that in FTC. In particular, TAMs were observed together with cancer cells in lymphatic invasive lesions, and the number of M2 cells inside and outside the lymphatic vessels showed a significant correlation. MMP-2 mRNA was expressed in nonneoplastic stromal cells as well as cancer cells, and double immunofluorescence staining confirmed M2 positivity. Consequently, this study reveals that M2 TAMs around lymphatic vessels within the tumour border of PTC may be associated with the lymphatic invasion of cancer cells. This study represents a step forward in elucidating the mechanism of lymphatic invasion.

A new risk factor indicator for papillary thyroid cancer based on immune infiltration.

Increasing evidence has indicated a close association between immune infiltration in cancer and clinical outcomes. However, related research in thyroid cancer is still deficient. Our research comprehensively investigated the immune infiltration of thyroid cancer. Data derived from TCGA and GEO databases were analyzed by the CIBERSORT, ESTIMATE, and EPIC algorithms. The CIBERSORT algorithm calculates the proportions of 22 types of immune cells. ESTIMATE algorithm calculates a stromal score to represent all stromal cells in cancer. The EPIC algorithm calculates the proportions of cancer-associated fibroblasts (CAFs) and endothelial cells (ECs), which are the main components of stromal cells. We analyzed the correlation of immune infiltration with clinical characteristics and outcomes of patients. We determined that the infiltration of CD8+ T cells improved the survival of thyroid cancer patients. Overexpression of immune checkpoints was closely related to the development of thyroid cancer. In general, stromal cells were associated with the progression of thyroid cancer. Interestingly, CAFs and ECs had opposite roles in this process. In addition, the BRAFV600E mutation was related to the upregulation of immune checkpoints and CAFs and the downregulation of CD8+ T cells and ECs. Finally, we constructed an immune risk score model to predict the prognosis and development of thyroid cancer. Our research demonstrated a comprehensive panorama of immune infiltration in thyroid cancer, which may provide potential value for immunotherapy.

OGDHL closely associates with tumor microenvironment and can serve as a prognostic biomarker for papillary thyroid cancer.

BACKGROUND: Papillary thyroid cancer (PTC) is the most common type of thyroid cancer. However, due to the lack of reliable prognostic biomarkers for PTC, overtreatment has been on the rise. Therefore, our research aims to identify new and promising prognostic biomarkers and provide fresh perspectives for clinical decision making. METHODS: The RNA-seq data and clinical data of PTC samples were obtained from The Cancer Genome Atlas data portal. GSE64912 and GSE83520 datasets were downloaded through the GEOquery R package. The difference in the expression of oxoglutarate dehydrogenase like (OGDHL) between PTC and normal tissues was explored by the Wilcoxon test. Kaplan-Meier (KM) and Cox regression analyses were used to further explore the prognostic value of OGDHL. The tumor microenvironments of PTC patients were explored based on ssGSEA and Tumor Immune Estimation Resource online database. Gene Set Enrichment Analysis (GSEA) was performed to explore the biological processes associated with OGDHL. RESULTS: The expression level of OGDHL in PTC was significantly altered compared to that in normal tissues (p < 0.05). Various biological processes associated with OGDHL were also explored through GSEA. KM analysis suggested that the low-OGDHL group had a better overall survival [OS, p = 3.49e-03, hazard ratio (HR) = 4.567]. The receiver operating characteristic curve also indicated the favorable prognostic potential of OGDHL. Moreover, OGDHL was proved to be an independent prognostic indicator in Cox analysis (p = 1.33e-02, HR = 0.152). In the analysis of the tumor microenvironment, the low-OGDHL group showed a lower immune score and stromal score, while tumor purity was higher. The expression of OGDHL was also closely correlated with the infiltration of immune cells. CONCLUSION: Our study elucidated the influence of OGDHL on the prognosis of PTC and demonstrated its potential as a novel biomarker, which would provide new insights into the prognosis monitoring and clinical decision making in PTC patients.

Predictors and a Prediction Model for Central Cervical Lymph Node Metastasis in Papillary Thyroid Carcinoma (cN0).

Objectives: To screen out the predictors of central cervical lymph node metastasis (CLNM) for papillary thyroid carcinoma (PTC) and establish a prediction model to guide the operation of PTC patients (cN0). Methods: Data from 296 PTC patients (cN0) who underwent thyroid operation at the Second Affiliated Hospital of Chongqing Medical University were collected and retrospectively analyzed. They were divided into two groups in accordance with central CLNM or not. Their information, including ultrasound (US) features, BRAFV600E status, and other characteristics of the two groups, was analyzed and compared using univariate and multivariate logistic regression analyses, and the independent predictors were selected to construct a nomogram. The calibration plot, C-index, and decision curve analysis were used to assess the prediction model's calibration, discrimination, and clinical usefulness. Results: A total of 37.8% (112/296) of PTC patients had central CLNM, and 62.2% (184/296) did not. The two groups were compared using a univariate logistic regression analysis, and there were no significant differences between the two groups in sex, aspect ratio, boundary, morphology, hypoechoic nodule, thyroid peroxidase antibody, or tumor location (P>0.05), and there were significant differences between age, tumor size, capsule contact, microcalcifications, blood flow signal, thyroglobulin antibodies (TgAb), and BRAF gene status (P<0.05). A multivariate logistic regression analysis was performed to further clarify the correlation of these indices. However, only tumor size (OR=2.814, 95% Cl=1.634~4.848, P<0.001), microcalcifications (OR=2.839, 95% Cl=1,684~4.787, P<0.001) and TgAb (OR=1.964, 95% Cl=1.039~3,711, P=0.038) were independent predictors of central CLNM and were incorporated and used to construct the prediction nomogram. The model had good discrimination with a C-index of 0.715. An ROC curve analysis was performed to evaluate the accuracy of this model. The decision curve analysis showed that the model was clinically useful when intervention was decided in the threshold range of 16% to 80%. Conclusion: In conclusion, three independent predictors of central CLNM, including tumor size (> 1.0 cm), US features (microcalcifications), and TgAb (positive), were screened out. A visualized nomogram model was established based on the three predictors in this study, which could be used as a basis of central cervical lymph node dissection (CLND) for PTC patients (cN0).

Inhibition of BRAF Sensitizes Thyroid Carcinoma to Immunotherapy by Enhancing tsMHCII-mediated Immune Recognition.

CONTEXT: Multiple mechanisms play roles in restricting the ability of T-cells to recognize and eliminate tumor cells. OBJECTIVE: To identify immune escape mechanisms involved in papillary thyroid carcinoma (PTC) to optimize immunotherapy. SETTING AND DESIGN: iTRAQ analysis was conducted to identify proteins differentially expressed in PTC samples with or without BRAFV600E mutation. Molecular mechanisms regulating tumor cell evasion were investigated by in vitro modulations of BRAF/MAPK and related pathways. The pathological significance of identified tumor-specific major histocompatibility complex class II (tsMHCII) molecules in mediating tumor cell immune escape and targeted immune therapy was further evaluated in a transgenic mouse model of spontaneous thyroid cancer. RESULTS: Proteomic analysis showed that tsMHCII level was significantly lower in BRAFV600E-associated PTCs and negatively correlated with BRAF mutation status. Constitutive activation of BRAF decreased tsMHCII surface expression on tumor cells, which inhibited activation of CD4+ T-cells and led to immune escape. Pathway analysis indicated that the transforming growth factor (TGF)-ß1/SMAD3-mediated repression of tsMHCII, which could be reversed by BRAF inhibition (BRAFi). Targeting this pathway with a combined therapy of BRAF inhibitor PLX4032 and anti-PD-1 antibody efficiently blocked tumor growth by increasing CD4+ T-cell infiltration in a transgenic PTC mouse model. CONCLUSIONS: Our results suggest that BRAFV600E mutation in PTC impairs the expression of tsMHCII through the TGF-ß1/SMAD3 pathway to enhance immune escape. Combined treatment with PLX4032 and anti-PD-1 antibody promotes recognition and elimination of PTC by the immune system in a pre-clinical mouse model, and therefore offers an effective therapeutic strategy for patients with advanced PTC.